function forth( X,Y,options)
q = forsythRegression(X,Y,str2double(options.p));
plot(X,polyval(q,X),options.color);
end

